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Tails of the Unexpected: The Interaction of an Isothermal 
Shell with a Cloud 
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ABSTRACT 

A new mechanism for the formation of cometary tails behind dense clouds or globules 
is discussed. Numerical hydrodynamical models show that when a dense shell of swcpt- 
up matter overruns a cloud, material in the shell is focusscd behind the cloud to form 
a tail. This mode of tail formation is completely distinct from other methods, which 
involve cither the removal of material from the cloud, or shadowing from a strong, 
nearby source of ionization. This mechanism is relevant to the cometary tails seen in 
planetary nebulae and to the interaction of superbubble shells with dense clouds. 
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1 INTRODUCTION 

An extensive literature of analytical and numerical inves- 
tigations of shock-cloud interactions now exists in which 
the effects of magnetic fields, radiative cooling, thermal 
cond uction and turbulence have all been considered (see , 
e.g., iPittard et all 120091 ; iPittard. Hartquist fe FaUel l20ld 
and references therein). However, in many astrophysi- 
cal objects the cloud size is comparable to or larger 
than the depth of the post-shock layer, so it is surpris- 
ing that to date almost all numerical investigations have 
been performed in the small-cloud limit where the post- 
shock flow has effectively an infinite depth. Notable excep- 
tions include investigations where the global flow is simu- 
lated, for instance, in the interaction of a supernova rem- 
nant (SNR) with a cloud (iTenorio-Tagle fc R6zvczkalll986l ; 



iRozvczka fc Tenorio-Tagld Il987l ; iLeao et al.ll2009l ), though 
such models suffer from low numerical resolution. Higher 
resolution was used in work on the t riggered collapse of a 
cloud by a thin shell (|Boss et al.ll201ol ). but the downstream 
flow was not examined. In this paper we re-examine the in- 
teraction of a cold, dense, isothermal shell with a spherical 
cloud, studying the effect of the shell thickness and Mach 
number. 



2 SHELL-CLOUD INTERACTION 

The shell is driven by a high pressure, low density, hot bub- 
ble. Pressure-driven shells occur in planetary nebulae (PNe) 
and around individual stars or groups of massive stars, and 
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have diameters ranging from ~ 0.1 pc to a few kpc. Since 
cooling breaks the scale-free nature of adiabatic simulations 
one is forced to choose a particular lengthscale. In this pa- 
per the interaction of a superbubble with a molecular cloud 
is simulated, though the results are also relevant to the 
smaller-scale interactions mentioned above. 

The structure of the superbubble, including the thick- 
ness of the shell and the density contrast between the shell 
and the hot interior, can be obtained from a global hydrody- 
namical simulation. However, the shell properties depend on 
assumptions about the structure of the ISM (e.g., whether 
there is a density gradient, the strength and orientation of 
magnetic fields, the number, size and density distributions 
of clouds, etc.) and the efficiency of mass transfer from the 
cool shell and from clouds within the hot bubble to the hot, 
rarefied gas. Therefore, we choose instead to specify instan- 
taneous shell and bubble properties (Mach number, surface 
density, bubble temperature) in order to focus on the key 
physics of the interaction. The shell is assumed to be planar 
(i.e. its radius is much larger than the cloud). 

Our investigation is based around a standard model 
with the following parameters. The inter-cloud ISM has 
an assumed density and temperature of p am b = 3.33 x 
1(T 25 gcm~ 3 and T amb = 8000 K. The density contrast 
of the spherical cloud x = 10 3 , and its core temperature is 
T c = 8 K. The cloud is thus in pressure equilibrium with 
its surroundings. To keep matters as simple as possible in 
this first investigation all material is assumed to behave as 
an ideal gas with a ratio of specific heats, 7 = 5/3, al- 
though in reality the cores are molecular. This simplification 
has little effect on the overall dynamics. The ISM pressure 
Pamb = 2667 cm" 3 K. The cloud has a mass M c = 300 M Q 
and a radius r c = 2pc. The density within the hot bubble 
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pbub = Pamb/100 ensures that this material does not appre- 
ciably cool during the simulations. We ignore thermal con- 
duction, magnetic fields and self-gravity. Though the latter 
is important for the parameters we have chosen, our focus in 
this paper is on the development of a tail behind the cloud, 
rather than star formation within it. 

The net heating/cooling rate per unit volume is param- 
eterized as e = nT — n 2 A, where n = p/mrr = 1.43 bh, 
and F and A are heating and cooling coefficients which are 
assumed to depend only on temperature. In the ISM, Y de- 
creases with increasing density as the starlight, soft X-ray, 
and cosmic ray flux are attenuated by the high column den- 
sities associated with dense clouds. Because the exact form 
of the attenuation depends on details which remain uncer- 
tain (for instance the size and abundance of PAHs), the 
heating rate at T < 10 4 K is similarly uncertain. In this 
work we assume that T = 10 _26 ergs _1 (independent of 
p or T). The low-temperature (T < 10 4 K) cooling was 
then adjusted to give 3 thermally stable phases at thermal 
pressures between 2000 — 6000 cm -3 K, as required by ob- 
servations. These stable phases, at temperatures ~ 10 K, 
~ 150 K and ~ 8500 K, correspond to the molecular, atomic 
and warm neutral/ionized phases, respectively. The cooling 
curve and phase diagram are shown in Fig. [T] The heating 
and cooling rates at low temperatures are higher than in re- 
ality, but since low temperature material is usually at high 
density (to maintain pressure balance) the timescale to reach 
equilibrium is short, and the simulations presented are not 
expected to be significantly affected by these simplifications. 

The pressure of the shell is high enough that the gas 
within it should be in the cold atomic or molecular phase. 
However, in all the simulations the column density of the 
shell is much less than that of the molecular cloud, so it is 
assumed that the heating is high enough to keep the shell 
material in the warm ionized phase. Therefore, the tempera- 
ture of material in the shell is prevented from cooling below 
8000 K, and the shock sweeping through the ISM behaves 
isothermally. In contrast, material ablated from the cloud is 
allowed to find its own equilibrium temperature, which may 
be higher or lower than T am b. 

W e adopt the density profile noted in IPittard et alj 
(|2009h with pi = 10 (i.e. a cloud with a reasonably sharp 
edge). Although the cloud is initially in pressure equilibrium 
with its surroundings, the adopted cooling curve means that 
the edge of the cloud is not in a stable thermal configura- 
tion. However, the shell hits the cloud before it has time to 
substantially adjust itself, and the initial transient is unim- 
portant to the following results. In contrast, the intercloud 
ISM is close to thermal equilibrium and its temperature re- 
mains almost constant throughout the simulation. 

A key parameter for the destruction of non-magnetic, 
non-conducting, non-gra vitating adiabatic clouds is the 
cloud-crushing timescale (I Klein. McKee fc Colellalll994r i. 

where v B and M a are the speed and adiabatic Mach num- 
ber of the shock through the intercloud medium having an 
adiabatic sound speed a. Although our study involves ra- 
diative shocks, the above characteristic timescale remains 
useful and is used in our analysis. All times are measured 
relative to when the shell is level with the cloud centre. 
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Figure 1. Top: The parameterized cooling curve adopted in this 
work. Bottom: The density-pressure phase diagram corresponding 
to the cooling curve shown and T = 10 -26 ergs -1 . 

3 NUMERICAL SIMULATIONS 

We use the MG code recently developed by Falle. It employs 
an exact Riemann solver for gasdynamics and piece- wise lin- 
ear interpola tion to ach ieve second order accuracy in space 
and time fcf. lFallelll99ll ). Adapative mesh refinement is han- 
dled on a cell-by-cell basis, and it is fully parallelized using 
MPI. Physics modules, including ones to handle self-gravity 
and magnetic fields, can be turned on or off as required. 
One such module is a k-e subgrid model, which we use in 
this work to calculate the properties of the turbulence which 
is generated in high Reynolds number astrophysical flows. 
Such models simulate the effect that the turbulent eddies 
have on the mean flow by increasing the transport coeffi- 
cients, particularly the viscosity, in regions where there is 
a lot of turbulence. Further deta ils including the full set of 
equations solved can be found in IPittard et all l|2009h . The 
2D calculations are computed on an rz axisymmetric grid, 
with a domain of ^ r ^ 24, —60 ^ z ^ 100 for all models. 
8 grid levels are used with 128 cells per cloud radius on the 
finest grid. A 3D calculation with a domain of —30 ^ x Si 30, 
0^y^3,0^;z^3 and 64 cells per cloud radius was also 
performed. The cloud is initially centered at the origin. 

Fig. [2] shows the time evolution of the interaction of 
a 2D model with an isothermal shell of thickness l s h = 8 pc 
and Mach number M a = 1.5. The pressure and density jump 
in the shell is 7M 2 = M 2 = 3.75 (where M a and Mi are re- 
spectively the adiabatic and isothermal Mach numbers of 
the shell). The ratio of the column density through the shell 
and the centre of the cloud, cr s h/cr c i = 7.5 x 10~ 3 . Hence the 
cloud is relatively unaffected by the initial passage of the 
shell, though the shocks driven into it by the jump in the 
external pressure start to compress it. This compression is 
nearly isotropic, due to the low Mach number of the inter- 
action and the high density contrast of the cloud. The cloud 
then exits through the back surface of the shell to reside in 
the low density interior of the hot bubble driving the shell. 

The most interesting aspect of the interaction is the 
formation of a tail behind the cloud (for earlier simu- 
lations which display similar but much thinner tails see 
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Figure 2. The time evolution of a 2D rz-axisymmetric shell-cloud interaction with M = 1.5, X = 10 3 an d an initial shell thickness of 
8pc (4r c ). The interaction proceeds left to right with i = 0.098, 0.162, 0.227, 0.291, 0.355, 0.549 and 0.806 t cc . Material in the shell is 
focussed into a tail as it passes over the cloud which is compressed due to the jump in pressure. The density maps span over 5 dex: the 
bubble interior has njr = 1.5 X 10 — 3 cm -3 (white), the ambient medium, shell and tail have njr = 0.15 — 0.6 cm -3 (light grey), and 
the unshocked cloud has njj = 150 cm -3 (dark grey). The highest speed attained by the flow is 23.8 kms" 1 . The unit of length is 2pc. 
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Figure 3. As Fig. [2]but showing the log of the pressure (black is high) with contour levels every 0.1 dex. The pressure gradients which 
exist around the cloud in the initial stages of the interaction focus material from the shell onto the symmetry axis. 



ITenorio- Tagle fc R6 Z yczkalll984bl ). The tail is mainly com- 
posed of material from the shell, with only small amounts 
(less than a few percent concentration) of material ablated 
or stripped from the cloud. The part of the shell adjacent 
to the cloud moves in the lateral direction onto the axis 
due to the pressure gradient which exists across its face as it 
sweeps over the cloud (see Fig. [3]). A large eddy forms which 



causes this material to lose speed relative to the rest of the 



the rear of the cloud (see also 


Nittmann. Falle & Gaskelll 


1 19821; ITenorio- Tagle & Rozvczka 


1984al). This material is 



subsequently compressed against the axis by the hot, sub- 
sonic flow which overtakes it. The pressure gradient across 
the face of the shell diminishes as the shell moves further 
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Figure 4. Same as Fig. [2] but showing cuts through the 2 = 
plane of a 3D simulation at t = 0.806 t cc - The shell is moving 
in the +x direction, and there is reflection symmetry about the 
y = plane. The panels show logarithmic values (black is high) 
of: a) the mass density and icy-velocity vectors; b) the turbulent 
energy density per unit mass; c) the temperature. 



creasing the thickness of the shell enhances the stripping of 
material from the cloud, which causes oscillations in the tail 
width (Figs.[6]3 and c). Interactions at higher Mach number 
enhance the growth of instabilities in the shear layer sur- 
rounding the tail (Fig. |6ji) . A model with a lower cloud den- 
sity contrast (x = 125 instead of 10 3 ) still produces a tail 
(Figs. [6^-g). Tails still form behind clouds with smoother 
density profiles and when the shell is curved (not shown). 



4 DISCUSSION 

The crushing of clouds by 
has been investi gated only 
the literature 



isothermal shells 

few times in 

Tenorio-Taele fc Rozvczkal Il984bl : 



iRozvczka fc Tenorio- Taeld 1 19871 : iLeao et all 120091 ) . While 
these works have demonstrated that a tail composed of 
shell material can form in an interaction with a cloud, the 
shell, and consequently the focussed tail, is always much 
thinner than the cloud radius. Furthermore, some of the 
tails are short-lived while others are soon dominated by 
ablated material. In contrast, we emphasize that the models 
presented here have long-lived tails of thickness comparable 
to the size of the cloud. 




Figure 5. Velocity along the core of the tail at t = 0.806 t cc for 
the interaction shown in Fig. [2] The cloud remains near z = 0.0. 
The subsonic flow in the interior of the bubble (from the left in 
this plot) decelerates as it approaches the cloud. 



downstream, and the focusing becomes more gradual. A 
3D simulation of this interaction reveals the same features 
(see Fig. [4|, indicating that this is a robust result which 
is not dependent on the assumed axisymmetry. The high 
shear around the cloud causes a turbulent boundary layer 
at the edge of the tail which grows with an opening angle 
of ~ 3 — 4° . The interior parts of the tail also contain some 
turbulence, though the central part of the cloud has none. 

The tails in these models exhibit a large length-to- width 
ratio which can reach nearly 50:1 at late times (t — 2.1 i cc ). 
Fig. [5] shows the velocity profile along the symmetry axis 
through the core of the tail at t = 0.806 t cc for the interac- 
tion shown in Fig. The acceleration is approximately con- 
stant along the tail, with the velocity reaching « 10 kms -1 
(M a = 0.8) at 2 = 10 r c . Due to the lack of material be- 
ing stripped off the cloud, the tail eventually dissipates as 
it thins and then detaches from the cloud (by t ~ 3t C c)- 
The axial velocity profile perpendicular to the tail shows a 
near constant speed at a given downstream position, which 
indicates efficient momentum transfer across the tail. 

We have performed a series of models designed to ex- 
plore parameter space to determine the conditions necessary 
for tail production (see Fig. |6j . Decreasing the thickness of 
the shell (i.e. er s h/(T c i) leads to a thinner tail (Fig. |6ji). In- 



4.1 Planetary Nebulae 

Cometary tail-like structures, formed behind dense molecu- 
lar clouds, are seen in many PNe. The best studied are those 
of the highly evolved PNe NGC 7293 (the Helix Nebula). The 
clouds are ionized on the parts of their surfaces exposed to 
the ionizing flux from the central star, and the tails point ra- 
diall y away from it. The ta ils often contain molecular mate- 
rial (|Matsuura et alj|2009l . e.g.,). Two different models have 
been proposed to explain the tails. In "shadow" models the 
tail forms due to the shielding of the direct ionizing radia- 
tion field of the ce ntral star (e.g.. lLopez-Martin et al1l200ll ; 
lO'Dell et ai]l2005l ). In contrast, "stream-source" models as- 
sume that the tail forms from material photoablated from 
the cloud (iDvson. Hartquist fc Bird 1 19931 ; iFalle et alj|2002l ; 
iPittard et al.ll2005l ; iDvson et alj|2006t ). 

The correct model is still dispute d (e.g., IDvson et al.l 
l2006t lO'Dell. Hennev fc Ferlandl 120071 ). However, observa- 
tions of the dynamics of the tails favour stream-source mod- 
els: i) there is no evidence for significant ionized ga s veloci- 
ties perpendicular to the tails ( Mcabur n et al ]| 19981). in con- 
trast to the shadow model of lLopez-Martin et al. r feoOlf ); 
the flow accelerates along the tails (by about 8—14 kms -1 , 
see iMeaburn fc BoumislBoTbT ). While our models are not 
specifically of the Helix tails, the velocity increase along the 
tail is similar to the observations, which are suggestive of a 
linear velocity gradient (Meaburn, private communication). 



4.2 Superbubbles, Starbursts and Superwinds 

The interaction of flows with clouds may be a key mech- 
anism for producing the broad emission wings to Ha 
line profiles seen in regions containing concentrations of 
massive stars including super star cluste rs and many giant 
HII regio ns, such as 30 Doradus (IChu fc Kennicuttl 
1994 iMelnick. Tenorio- Tagle fc Terlevichl Il999h . 
Westmocmette et al.l ( 2007a! bl) concluded that the broad 




Figure 6. (a) M = 1.5, l ah = 2pc (r c ), t = 0.806 t cc . (b) and (c): M = 1.5, l sh = 32 pc, t = 0.806 and 1.192 t cc . (d) M = 3.0, l ah = 8pc, 
t = 0.519 tcc (e)-(g): X = 12 5, M = 1.5, Z sh = 2pc (r c ), i = 0.642, 1.006, and 1.552 t cc . Note that the cloud is thermally unstable in this 
model, and collapses to higher densities and much colder temperatures than its initial conditions. The density grey-scale is consistent in 
all panels, but the velocity vectors in panel d) are twice as short for a given speed compared to the other panels. 



component arises in a turbulent boundary layer at the in- 
terface between hot gas flowing past cold gas stripped from 
clouds. However, subsequent modelling has i ndicated that 
unrea listically high flow speeds are needed l|Binette et al.l 
l2009h . Instead, we believe that the broad emission wings 
reflect the acceleration of material along the tail and not 
just turbulent motions within the mixing layer. 

Beautiful filamentar y structures are a ls o seen in star- 
burst superwinds (e.g., ICecil et al.l 1200 ll ; lOhvama et alj 
120021 ). While these may be materi al stripped from d enser 
clouds, as seen in simulations (e.g.. ICooper et al]|2008l ). the 
results in this paper indicate that some of the filaments may 
in fact have been formed directly out of an overrunning shell. 



5 CONCLUSIONS 

We have demonstrated a new mechanism for the formation 
of tails behind dense clouds which involves the removal and 
trailing of material from an overrunning isothermal shell. 
The mechanism appears robust to a range of shell thick- 
nesses and Mach numbers, and the cloud density contrast, 
though these parameters influence the tail's properties. 
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